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One- and two-photon transitions in the hydrogen atom are analytically evaluated in the absence and in the 
presence of an external electric field. The emission probabilities are different for the hydrogen (H) and anti- 
hydrogen (H) atoms due to the existence of contributions, linear in electric field. The magnitude of these 
contributions is evaluated within the nonrelativistic limit. The Coulomb Green function method is applied. 
Different nonrelativistic "forms" for the decay probabilities in combination with different gauge choices are 
considered. The three-photon ElElEl 2p-ls transition probability is also evaluated and possible applications of 
the results are discussed. 
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INTRODUCTION 

Recent experimental success in the production of anti-hydrogen atoms opens reahstic chances for a comparison of the 
atomic spectra for hydrogen (H) and anti-hydrogen (H) atoms. One of the purposes for this comparison is the search for the 
CPT-violating effects fl- Thepossibility of the CPT-tests is connected with the modern extra- accurate resonance frequency 
measurements in hydrogen ['41,1511 • In l@l-fl it was shown that a specific difference in the H and H atomic spectra arises even in 
the absence of the CPT-violation, if an external electric field is present. In principle, a difference arises also for the frequency 
measurements if nonresonant (NR) coiTections are taken into account. Actually, the NR coiTections define the natural limit 
up to which the frequency measurements have sense: beyond this limit the spectral line profile cannot be defined by the two 
parameters: resonance frequency and width. NR connections are not just corrections to the level energy, but rather coiTections to 
the spectral line profile. These corrections for the atoms placed in an external electric field contain terms, linear in the electric 
field. Therefore NR corrections appear to be different for H and H atoms. However, this difference is very small and hence the 
existence of the electric stray field should not become a serious problem in performing experiments searching for CPT-violating 
effects. Another difference in spectroscopic properties of H and H atoms in an external electric field is the difference in the 
transition probabilities, also caused by the terms, linear in electric field. This effect is much larger than the difference in the NR 
coiTections and was not yet discussed in detail. 

In the present paper we will investigate specially the difference in the one-photon 2s- Is emission probabilities in H and H 
atoms in an external electric field and will find the optimal conditions for observing this effect. This investigation should confirm 
our understanding of fundamental symmetries in physics. Thus, going over from H to H and changing simultaneously the sign 
of the electric field should not change the atomic spectra provided that CPT symmetry is conserved. 

The one-photon decay of the 2s state of the H atom in an external electric field was utilized in the experiments 10, Btl for the 
registration of the 2s- Is two-photon resonance absorption. The difference of the one-photon 2s- Is decay rate for H and H atoms 
in an external electric field can be observable in the similar experiments with H atoms. 

The probabilities for the spontaneous two-photon decay in hydrogen atoms and hydrogen-like ions are under investigation 
since the theoretical formalism has been introduced by Goppert-Mayer ISJ and the first evaluation for the two-photon ElEl 
transition 2s 27(E1) + Is has been presented by Breit and Teller lllOll . A highly accurate calculation of the ElEl - transition 
probability has been performed by Klarsfeld flTl. Recently Jentschura fl^l performed a complete evaluation of the radiative 
corrections and presented more accurate value of the ElEl two-photon decay probability. The double- and triple-photon decays 



of metastable "'Pq atomic state were considered in 11311 . The present paper is devoted also to evaluation of the probabilities for 



two-photon decays 2p 7(E1) + 7(M1) + Is and 2p 7(E1) + 7(E2) + Is. Evaluations of these two transitions have 



been first accomplished in II141I15I1 for hydrogen-like systems with nuclear charge numbers Z within the range 1 < Z < 100 by 
pure numerical methods. Here we present analytic calculations in the nonrelativistic limit and compare them with corresponding 
numerical results. For performing the summations over intermediate states (i.e. over the complete set of solutions of the 
Schrodinger equation describing electrons in the Coulomb field of the nucleus) we employ the Coulomb Green function lfl6ll . 
The Green function method has been first applied for deriving the general expression for the two-photon decay probability in 



H atom and H-like ions in lllTl [Is]. An alternative approach applicable for arbitrary states based on Schwinger's analytical 
representation of the Coulomb Green function has been presented in 11 9, ,20.1 . 

As it has been indicated in 121], the nonrelativistic behavior of ElMl transitions as a function of Z with the neglect of the 
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interelectron interaction should be Vl^^^^i = (8/97r)(aZ)^^/100 for the helium-Hke systems. This very small value arises 
due to the cancellation of contributions of the leading terms 2pi/2 and 2p3/2 in the summation over intermediate np-states. 
However, we should note that this result yields only a minor contribution for small nuclear charges Z when it is evaluated within 
the "velocity" gauge 12 111. In this case a major contribution arises from the negative-energy intermediate states and scales like 
{aZY in atomic units Hj], iflsll . 

The two-photon 2p]^^2 ~^ 7(E1) + 7(E2) + Is transition rate for hydrogen and hydrogen-like ions is proportional to (ctZ)^ 
in atomic units. Thus the two-photon transitions represent higher-order corrections to the life time of the 2pi /^.-level when 
compared to the lowest order {aZ)'^ \n{aZ) (in relativistic unit) radiative corrections derived in references 12211 . 1230 . A direct 
observation of the influence of the two photon 2pi/2 ~^ 7(E1) + 7(E2) + Is transitions in the H atom does not look feasible 
due to the huge background arising from the one-photon transition 2pi^2 ^ 7(E1) + Is. However, two-photon decays of the 
2pi/2-level could be observed in coincidence experiments. 

In this paper we reevaluate the two-photon decay rates of the 2p-state in hydrogen, deriving the E1E2 and ElMl contributions 
to the two-photon emission processes by means of analytical methods. These calculations are performed within different gauges 



and employing different forms for the expression of the transition probability (see 112411 . II25I1 V 

We also perform the calculation of two-photon transition probabilities from the 2s and 2pi/2 hydrogenic levels in the presence 
of an external electric field. The mixing of the 2s and 2pi/2 levels results in additional terms, linear in electric field. The 
evaluation is performed within the nonrelativistic limit by the Coulomb Green function method lIlSll in the "length" gauge. 



These terms linear in the electric field lead to the difference of the corresponding probabilities, see Il26ll - ll28ll . for the hydrogen 



and anti-hydrogen atoms. The magnitude of these contributions is evaluated and the possibility of the observation of this effect 
is discussed. It is important to stress, that these linear field corrections even in a very small stray electric field can be larger than 



the radiative correction to the decay of 2s level, evaluated in Ill2ll . 



Finally, we evaluate a three-photon emission probability for the hydrogen atom. This question was never debated in literature. 

The recent success in observation of the cosmic microwave background temperature and polarization anisotropy draw attention 
to the details of the cosmological hydrogen recombination history. This, in turn, required the accurate knowledge of the two- 
photon decay processes in hydrogen (see |2^ for details and references). The main interest concentrates around the 2s- Is 



ElEl two-photon transition and the two-photon decays for the ns, nd excited states. 

The smallness of the numerical coefficients in E1E2 and ElMl two-photon expressions excludes any significant role of these 
transitions for astrophysics. The same can be stated for the 3-photon 3E1 2pi/2~ls decay. However the values for the 2 p i /o — Is 
ElMl and E1E2 transition rates appeared to be important for the determination of the theoretical accuracy limits 831113211 for 
frequency measurements by means of optical resonance experiments with hydrogen Hi. 

The following notations will be used throughout this paper. The 2pi/2 level is labelled as 2p because the 2p3/2 level will not 
be considered. Vectors with 3 components are in bold. Angular part of any vector r is defined by the unit vector rir directed the 



same as r. The Clebsh-Gordon coefficients Cl^_^ ^^^^ are introduced according to ll33ll . Other notations are the following: 



A = jAi?2 + -r^p, (1) 



^(17) ^ ^wi':^/wil-^\ (2) 



^(2.) ^ ^w!,T^/W^l'\ (3) 

Here AE^l is the Lamb shift, r2p is the width of the 2p state, VFa^^^^'^^^'' and W^'^''"'^'^'^ are the one-photon and two-photon 

decay probabilities for the corresponding 2s and 2p levels. The value w!^^ = w!^^^"^^ + VK2p ^^^^ includes the probabilities 

for 2p 7(E1) + 7(E2) + Is and 2p 7(E1) + 7(M1) + Is decays. 

2S DECAY RATE FOR HYDROGEN AND ANTI-HYDROGEN ATOMS IN EXTERNAL ELECTRIC FIELDS 

In the absence of an electric field, the Ml transition 2s —> Is + 7 - strongly forbidden in the nonrelativistic limit - was first 
evaluated by Breit and Teller [10]. The value, obtained in ifioll was later improved by Drake |34] (see also l,35ll '). In |34, 35] 



relativistic corrections to the Schrodinger wave functions and to the Ml photon emission operator were taken into account. 
Neglecting these corrections the 2s ^ Is + 7 transition probability turns to be zero due to the orthogonality of the radial wave 
functions. Accurate fully relativistic calculations for the 2s ^ Is + 7 transition for the hydrogen-like ions with arbitrary nuclear 
charge Z values were performed by Johnson ll36ll . 
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Within QED theory the S-matrix element for the one-photon emission process reads 1137 1 



Sa'A = ^ / dx^A'ix)e:re-''-''''-^'^Mx), (4) 



2uj 

where k is the wave vector of the photon, oj = |k| is the photon frequency, x — (r, t) is the 4-vector of space-time coordinates, 
is the 4-vector of the photon polarization, tpA (x) is the Dirac wave function for the bound electron in an atom, 7^ are the 
Dirac matrices and A', A correspond to the relevant quantum numbers characterizing the final and initial states of the electron 
in an atom. In Eq. (IDl e is the electron charge and the relativistic units h ^ c = 1 are used throughout this section. 
We integrate over the time variable t in Eq. ^ and employ the relation 

Sa'A - -27:iS{Ea' -Ea~ lo)Ua'a (5) 
for the transition amplitude Ua'A- In what follows, we define also the amplitude Ua'aO^, e) as 



V2 



TTLU 



Ua'a{^, e) = ^ Ua'a - (A'|eae-*''"-|^>- (6) 

e 

Here e is the transverse polarization vector and a are the Dirac matrices. 

Transition probability integrated over the photon emission directions and summed over the polarizations is 

wi]-'2 = '^j:J'i^-\UA'A{k,e)\\ (7) 

u!aa' — Ea — Ea', Ea', Ea are the Dirac eigenvalues for the atomic electron, is the matrix element with the Dirac 

4-component wave functions. Summation over e impUes the transversality condition ek = 0. 
Within the Pauli approximation Eq. ^ reduces to llssll : 

(8) 



with the transition amplitude in the Pauli approximation defined via the matrix element of the corresponding emission operator 

UIa{K e) = ((ep + ie[k x sDe"*"-)^,^ . (9) 

The relation Eq. (|9]l involves the nonrelativistic electron momentum operator p = — iV, the electron spin operator s ~ {a 
are the Pauli matrices) and {■■■)a'a = denotes the matrix element with the Schrodinger wave functions. The first 

term of the integrand in Eq. ^ describes the El photon emission, that in case of A' — Is, A — 2s is forbidden by parity; The 
second term corresponds to the Ml photon emission. Due to the orthogonality of the radial wave functions il'2s and ipis the Ml 
transition probability as described by Eq. ([8]) is nonzero only due to the factor e~^^'' and due to the relativistic corrections to 
the Schrodinger wave functions. For the atomic electron r ~ oq = 1/maZ (ao is the Bohr's radius, m is the electron mass, 
a « 1/137 is the fine structure constant), uj « m{aZ)'^ and kr « a. Thus, one can restrict the consideration to the lowest 
nonvanishing term of the multipole expansion of e^''"', which turns to be of the same order of magnitude as the relativistic 
corrections to the wave functions. 

The zero-order wave functions which should be employed for the evaluation of the matrix elements in Eq. ([8]) look hke 

ijAir) = V«,7™(r) = Rni{r)Yi^, (nr)x™., (10) 

where the standard set of one-electron quantum numbers is introduced: Principal quantum number n, total electron angular 
momentum j and its projection m, orbital electron angular momentum I and its projection nii, spin projection m^. The function 
Xm, is the Pauli spinor For the 2s and Is electron states I — and the transition amplitude in Eq. dH) reduces to I35i1 

Ulm^ism'^ = ~ik^{ms\e[k x a]\m'^)R2sis, (11) 

OQ 

i?2.is = J R2s{r)Ri,{ry dr, (12) 
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where |ms) are the spin wave functions, i?2s(?'), Ris{t) are the nonrelativistic radial wave functions. Inserting Eq. ( fTTT i into 
Eq. we have to sum over and average over nis- For summation over e the standard formula 

^e^ek ^ Sik - n,nk (13) 

e 

is used, where rii and Uk are the Cartesian components of the vector nk- After integration over the photon emission directions 
nk one finds 

W^isL^ = ^«"^'("^)''^2sl., (14) 

respectively. The evaluation of the radial integral in Eq. ( fT4l i results 

W^sll = ^"^«("^)'" = 2.8 • 10-6 s-i . (15) 

The decay rate of the 2s state of a hydrogen atom in the presence of an external electric field was studied in ll26ll -ll2 8| . The 
external electric field mixes the states 2s and 2p. A 100% mixing occurs in a field with the strength — 475 V/cm f38f|. In 
what follows we will consider weaker fields D < Dc such that admixtures of all other states, apart from 2p can be neglected. 
The ground state Is will be assumed unaffected by the field. We will denote the mixed states 2s and 2p as 2s and 2p, respectively. 
For the wave function of the 2s state we can write Il27il 

\2sm,) = \2ams) + ?7^(2prn;,'|eDr|2sms)|2pm;,'), (16) 

m'J 

where i] = (AE'l + ir2p/2)-^, D is the electric field vector. 

The transition amplitude Eq. (fTTT l in an external electric field looks like 

^Tsm.ism'S^^^) = f^2';™.isK(k.e) +^E(2s™s|eDr|2pml')C/2V.i,™, (k,e), (17) 

m'J 

where 

Uipm'j isr<(k,e) = iTOt^2sis(2pm'/|er|lsm',). (18) 
Direct evaluation of the integral in Eq. dTsl l results in 

Ulra'j is™i (k, e) - 3imco2s is E(-l)'e«^K< ' (l'^) 

where are the spherical components of the vector e and q = —q. Similarly, 

(2sm,|eDr|2pmi') = 3e E(-l)*'i?,Cfjr|^, . (20) 
g' 

Further evaluation requires the insertion of the amplitude Eq. ( fTTI ) in Eq. (|7]) and summation over e. For this purpose the 
formula Eq. ( fT3] ) should be used and then the representation of the scalar and vector products in spherical components should be 
employed. The final result is (here we do not integrate over the directions nj.): 



Tap e^D^ 
1 + eDnk- ^ 



i,(l7)A2 (w(l7)A)2 



dnk. (21) 



Although formula Eq. ( 1211 1 was obtained earlier in [261 and 12711 . 12811 . We present here our way for its derivation. 

The formal T-noninvariance of the factor Dnk in Eq. (1211 1 (nk and D are T-odd and T-even vectors, respectively) is compen- 
sated by the dependence on r2p; This is the imitation of T-noninvariance in unstable systems, as predicted by Zeldovich 13911 . 
Rewriting Eq. ( 1211 1 into the form ll26ll 

dW^^ll = T^o [1 T /3(i?)nDnk] dnk, (22) 
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where 



8^ 2sl« (y;(l7)A)2 



Wo = —W^p,i { 1 + ^^T^vTTT ) , (23) 



we find the maximum value (3{D) at ll26ll 



The value /J^ax = /3(£'max) is equal to 



/3(D) = '^'^^2pu;(^^) 



Anax = ~ 0.3 • 10-4 V/cm. (25) 

e 



/3max- 2A 20' ^^^^ 



The (— ) and (+) signs in Eq. ( |22] | correspond to the H and H atoms, respectively. 

The relative difference for the decay rates in H and H atoms at the maximum value Umax equals to: 

^<I>) _ ^^£is(H) _ Tyo(Anax)2;3(Jmax)nDni. 1 



In the presence of such a very weak field given by Eq. (1251 1. this difference is close to about 20% and probably can be observed 
in experiments of the type reported in |4, 5]. 

We would also note that if one integrates in Eq. ( |22] | over photon emission directions, the term hnear with respect to the field 
vanishes. But a quadratic term included in Wq exists and represents a correction to the W2s^7s transition probability, i.e. to the 
lifetime of the 2s level. This correction term reaches the magnitude of the radiative correction obtained in 1,121 in the field of the 
strength 



1 8^ sr,^ ^ 87: sr, 

^'-~Ti\hr^ — OVT V^"r — Anax, (28) 

1^1 V ^ ^2pTs V 3 r^s 

where the correction ST2s/^2s was derived by Jentschura (see Eq. (36) in [12]). Though the correction in flit] is obtained for the 
process of the two-photon decay of 2s level, it also represents a correction to the lifetime of the 2s level. For the hydrogen atom 
~ 1.4 • IQ-^ V/cm. Since the experiments deal with differential cross sections, we compare also the linear term with the 
radiative correction. The linear term /3(-D) reaches the magnitude of (5r2s/r2s at the field strength 7.5 • 10^^^ V/cm. It should 
be difficult to eliminate spurious fields of such magnitude in real experiments and, therefore, the comparison of the theoretical 
results in il2ll with experimental ones requires some caution. 



TRANSITION PROBABILITIES IN DIFFERENT FORMS AND GAUGES 

In this section different gauges in combination with different "forms" for the one-photon transition probability are described. 
Atomic units 7i = e = to = 1 will be used throughout this section. 

The transition probability for the emission of a photon with definite angular momentum and parity can be described in the 
first order of QED perturbation theory within arbitrary gauge as 



Wa-,A'{uj) = [K^'I KA^fc,(r)) + <i>^fe, (r) \A)f 



kq 

+ \{A'\a^-,A^kq {r)\A)f 



(29) 



where c-^ujkq and m-^ujkq denote the electric and magnetic vector potentials and ^ujkq corresponds to the scalar potential; oj 
is the photon frequency, k, q are the total angular momentum of the emitted photon and its projection. The bra-kets \A) and 
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{A'\ are stationary Dirac states (wave functions) with energies Ea and Ea' and a are the Dirac matrices. In the momentum 
representation these potentials take the form 



4^2^3/2 



3/2 



-(5 U 



(fc) = 



(30) 
(31) 
(32) 



Here k denotes the variable in the momentum representation. Functions o Y^.^ and mY^:^ are the vector spherical harmonics 
of electric and magnetic type, respectively, Y^q is the ordinary spherical harmonic, c is the speed of light and K denotes a 
gauge-dependent constant. 

The spherical components of the transversal electric c^'^^' and the longitudinal i^^^'*') parts of the electromagnetic vector 
potential (the superscript (lA) defines the rank of a spherical tensor and labels the components) are 



.A 



(lA) 

ujkq 



7rc(2fc + 1) 



+ v/(fc + l)(2fc- 1) 



1 fc- 1 fc 
X q — X q 



(-1) 



fc+g-A 



(33) 



Lukq 



OJ 



7rc(2fc + 1) 



+ V^(2fc-1) 



V(fc + l)(2fc + 3) 
A g-A l)9k-i{^r)C 



1 fc + 1 fc 
X q — X q 



(fc-l) 
g+A 



-k-l ^_^^k+q+X+l 



(34) 



where C 



(fc) 



47r 



'-9 Y 2/:+ 

magnetic vector potential look as 



(k) 

jYl^' and usual notations for 3j-symbols are employed. The spherical components of the transverse 



fc fc 

—A -q + X -q 



C 



(k) 

-g+A' 



while the spherical components of the scalar potential are given by the following expression 



^^kq = ^-\-l)''+'2J-g,iLur)Y^ 



(35) 



(36) 



The radial functions 5fc(wr) ai-e related to Bessel functions Jfj,{z) via gk{z) = (27r)'^/^-^ Jj,_|_ i (z). 

Usually two gauges are used: The so-called Coulomb gauge that corresponds to the vanishing longitudinal part of the vector 
potential and the scalar potential (i.e. V • A = V • o A = and <i> = 0). This gauge is characterized by the choice of the gauge 

parameter K = 0. Another convenient gauge is defined by the following value of the parameter K — ^\J~^^- Within this 

gauge, as it can be seen from Eqs. (1331 ) and ( |34| |. the terms containing spherical functions ci'^J of the order fc — 1, vanish in the 
expression for the transition probability (|29] l. 

After some manipulations the expression for the probabihty of emission of an electric photon with the angular momentum fc 
can be cast into the form 



W 



Ek 



2(fc + l)w 
fc(2fc + l)c 



E 

q— — k 



{A'\ 



/(fc) 



K 



k 



O 



(fc) 



\A) 



(37) 
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Here 



.0 



+ v/fc(2/c-l).gfe_i(cor) 



(fe) 



(fc) 



(fc) 



S.0 



(fe) 



v/(fc + l)(2fc + 3).gfc+iK) 



^(fe+i) X ad) 



(fc) 

-9 



X a 



(1) 



(38) 

(39) 
(40) 



and fa^i) X 

^(s2)j(«) represents the tensor product of two irreducible spherical tensors of rank si and S2 coupled to a spherical 
tensor of rank s with components q. 



Using the following integral relation for the Dirac wave functions ||37|] 



where x is an arbitrary function, one can establish another form for the Efc-transition probability (see [12511 ): 



W 



Efc 

A~>A' 



2(fc + l)a;3 
k{2k + l)c^ 



E 

q—~k 



(A'|eOL^)+X- 



uj\l k + 1 



(41) 



(42) 



where 



(fe) 



-9k i^^r) 



k + 1 

+ v/(fc + l)(2fc + 3) [C7(^-1) X ad) 



(fe)- 



C('=-i) X ad) 



(fc) 



(43) 



Thus, we have two different (equivalent) forms for the Efc-transition probabilities (Eqs. (|37] | and (021)) together with an arbitrary 
choice for the gauge constant K at our disposal. Analogous expressions for the emission probability of a photon, characterized 
by its energy and polarization, are provided in ll40ll . 

The aim of the present investigation concerns the derivation of the nonrelativistic limit of the E1M1-, ElEl- and E1E2- 
transition probabilities in different gauges and forms. Deriving the nonrelativistic limit of Eqs. ( |37] ) and (|42] |. implies two 
distinct nonrelativistic forms for the one-photon transition probability with arbitrary gauge constant K (see MM)' 



W 



Efe _ 2(fc + l)(2fc + l)cj2fe-i A 
^ fc[(2fc + l)!!]2c2fe+i 



k 



- -Q^'ii 1^) 



(44) 



and 



p^Efc 



A^A' 



2(fc + l)(2fc+l) 
fc[(2fc + 1)!!]2 



q— — k 



{A'\ 



Jk 



-Q 



/(fe) 



Q 



(fe) 



1^) 



Here \A) and {A'\ are nonrelativistic Schrodinger states (wave functions) together with operators 



Q 



(fe) 



knik) 



(45) 



(46) 



-fc-i / ,„Mk) ^ 

-q 



Ci^) X Ld) 



(fc) 



(47) 
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where i^^-' is the orbital angular momentum of the atomic electron. Choosing K = 0, we find that the operator in Eq. (|44] | 
corresponds to the nonrelativistic transition operator in the "velocity" form, while for K = —\ -^77^ it is related to the transition 



fc 

operator in the "length" form. However, the correspondence of a certain gauge choice to a particular type of nonrelativistic 
transition operators is not unique. In view of Eq. ( l45T l we can conclude, that within the nonrelativistic limit the expression 

Eq. ( l42b with K = Q converts the transition operator into the "length" form and with K = into the "velocity" form, 

respectively. 

APPLICATION OF THE COULOMB GREEN FUNCTION 

In order to evaluate the transition probabilities for the processes 2p — > 27 + Is and 2s ^ 27 + Is in the hydrogen atom the 
nonrelativistic Coulomb Green function is employed. The summations over the entire spectrum of the Schrodinger equation arise 
usually when perturbation theory is applied. The Green function approach allows one to express the intermediate summations 
in a closed analytic form. This is very useful for the analysis and for tests of numerical calculations. 

The Green function for the Schrodinger equation with the Hamiltonian H is defined by the solution of the equation 



- EjG{E;r,r') =S{r-r') (48) 
and can be always represented in terms of a spectral decomposition 

GiE,ry)^j:^^^^^. (49) 



N 



In Eq. ( |49] l the sum runs over the entire spectrum of the Hamiltonian (bound and continuous spectrum). The set of quantum 
numbers TV may be specified as usual by the principal quantum number n, orbital angular momentum number I and projection 
TO. In view of the spherical symmetry it is sufficient to derive a closed expression for the radial part gi{E; r, r') of the Green 
function defined by the partial wave decomposition 

G{E-v,v')^Y.^9i{E;r,r')YCMyirnM. (50) 

Im 

In the particular case of an external Coulomb potential the Green function of Eq. ( l48T l is called Coulomb Green Function (CGF). 
With the use of the expansion Eq. dSOl l the radial integrals occurring in Eqs. (l44l i and (|45] | for the transition probabilities can be 
calculated explicitly (see Ref. |42] for details). 

For the radial part of the Coulomb Green function it is convenient to employ the Sturmian expansion [18], which is defined in 
the entire complex energy plane via 



^ 00 
— gi{E,r,r')= ^ 



n^— 



(51) 



where <I>„^j (r) denote the Sturmian functions. The Sturmian expansion of the CGF can be represented in an alternative form by 
introducing radial functions 



which are related to Whittaker functions A/„ where = Z/ \/—2E. For integer values = n these functions coincide 

with the normalized, radial hydrogenic wave functions 

^nAr) = MRni(-\. (53) 



Z V V 



Substitution of Eq. dSSl l into dSTT i yields 

1 jy^ °° 



n 



Rnl ( —\ Rnl (— I . (54) 



n — V \ V \ V 
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Within this paper we apply the Green function method for the evaluation of the two-photon decay probability in the hydrogen 
atom. 

In 14311 the two-photon transition process 2s 7(E1) + 7(E1) + Is has been considered. The probability for the two-photon 
decay A 7(E1) + 7(E1) + A' with photon frequencies uji and uj2 within the nonrelativistic limit and dipole approximation 
is: 

dW^!^},{co,) = f (y)' E \iA'\rY,MAn^)GiEA-co,;vy)r%\,Jn^,)\A) 

^ ^ ' Ml Ah 

+ [A'lrYiMi (nr) G{Ea - u;2;r,r')r%\,^ (n,,) \A)\^ (cc-i^a)' dw2 ■ (55) 



The energy conservation law implies uji = Ea — Ea' — UJ2- After the evaluation of angular matrix elements in Eq. (1551 1 the 
remaining radial integrals have the form 



OO OO OO 



dr'drdx {r'y'+^ry+i exp ^--(/3V' + /3r) + (r + r')cosh(a;)^ 





X [coth y-jj hi+i I — ^sinh(a;) 1 , (56) 

where n and n' are the principal quantum numbers of the initial and final states, respectively, together with the parameters 
P ~ v/n, P' — v/n' and v = \/—2 (Eni — lo). These integrals can be evaluated analytically after inserting the series expansion 
for the modified Bessel functions l2i+i - The integration over x should be done at the end. Assuming emission of two El photons, 
one can write the total probability for such a two-photon decay as 



U)0 

WT^A' = \ I dWT_}A'{^) (57) 



with = Ea — Ea' ■ For the process 2s 27(E1) + Is in ll43ll the result of the evaluations was reported as 

Wi^^l = 8.226(aZ)6 s~^ (58) 

with an accuracy of about 1%. 

However, in further calculations it is convenient to utilize the other representation of the radial Coulomb Green function in 
terms of an expansion over Laguerre polynomials LI 8.1 (employed also in Ref. 1.43,1 ) 

= '^"H'^JS (2^ + l + n)!(n + / + l-.) - ^''^ 

The series Eq. (|59] l converges absolutely as n^'^l'^ for arguments r, r' > and Im(i^) = [ 1 8] . The angular momentum quantum 
number I ^ 1 for the intermediate states is fixed after the angular integration. Inserting the expansion ( |59] | for I = 1 into the 
expression dSSl l yields 

dWi^\{u;)=(iy^[h{.)+h{.')fa'dco, (60) 



CXJ 

\2/ ^-^ (m + 6)\(m + 2 — v) I 

m—O Q 

OO 

X J dU^e-'C^) (l - Ji) LUt). (61) 





These integrals can be evaluated analytically. After inserting this result into Eq. ( |57] | the integration over ui has to be performed 
in order to obtain the total transition probability for the ElEl decay of the 2s state in the hydrogen atom. This is achieved 
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numerically with the aid of the computer-algebra code MATHEMATICA. The final result is 

W^liii dW^^^^iuj) = 0.00131823 {aZf a.u. = 8.22932 s'^ {Z = 1) (62) 



with LUf) = i?2s ~ ^'is- In Eq- (l62] i we indicated the Z-dependence of the W^\f ^ transition probability. The numerical value 
(|62] | coincides with most precise result I,1Z1 up to 5 digits. This will serve us as accuracy estimate of our approach in following 
calculations. 



E1E2 DECAY PROBABILITY FOR THE 2P STATE 



In this section we consider the E1E2 decay of the 2p state in the hydrogen atom. Again the set of quantum numbers nlmi 
is employed as far as the total angular momentum j is not important in this calculation performed within the nonrelativistic 



approach. Nevertheless, in order to compare our results with those obtained from the relativistic evaluation (see Refs. 1114 1 



il5\\ ). we shall perform the calculation within two different gauges according to Eqs. (I44l l-(|47]|. This will also elucidate the 
potential influence of relativistic effects associated with the contribution of the negative-energy Dirac spectrum. 

As a test for the method the gauge constant K ~ — \f^^ chosen in the expression Eq. (l44l i for the transition probability, 
which corresponds to the nonrelativistic "length" form as mentioned above. This would be equivalent to the choice K — 
together with the form Eq. ( |45] ). Inspection of Eqs. ( |44| ) and ( l46l l reveals that the emission of electric photons (E/c) is described 
by the potentials 

fc+ 1 2uj^+^ 
~ {2k + l) 



Accordingly, the two electric photon decay rate of the atomic state A can be written as 

2 

S {uj + Lu' ~ Ea + Ea') dLudLu'.(64) 



^ {A'\V^''\N){N\V^'''\A) ^ (A'|yEfc'|7v)(7v|yEfe|^) 
— f: A -4-1,! ^ 



^ En-Ea+uj ^ En~Ea+uj' 



Here the labels A, A' and N abbreviate the set of nonrelativistic quantum numbers (principal quantum number n, orbital mo- 
mentum / and projection to/) for indicating the state of the atomic electron as the initial (A), intermediate (N) and final (A'). 
The photons will be characterized by the angular momentum and its projection (kq) as well as by the type of the photon (electric 
or magnetic). The Eq. ( |64| ) also implies the summation over degenerate substates of the final atomic state A' and the average 
over the degenerate substates of the initial atomic state A as well as summations over the angular momentum projections of both 
emitted photons. The frequencies of the two photons uj and uj' are related by the energy conservation law uj' = ujq — uj, where 
u;o = Ea - Ea- ■ 

Employing the eigenmode decomposition of the Coulomb Green function Eq. dSOl l the probabiUty of the two-photon decay 
process takes the form 



dW 



EfeEfc' 



27r 



T ^ 



21a , ^ , 



J2j J ^i"i'^i'2i?n^,;^,(ri)y/*,„^,(nrj 

Inii 



+ I '^^i'^^^ i?„„;,,(ri)y;*,„^,(nrjF^='='(ri)5;(^.';ri,r2)i1™,(nrJ 

Imi 

(65) 

where V^'' (r) is the potential (l63T l written in the gauge K = — -^ii and compatible with the form Eq. (l44l l. together with 

parameters v = Z/ \J —2(Ea ~ w), v' = Zj \J —2{Ea — t^') and frequency lo' = Ea — Ea' — lo. 

Specifying Eq. ( |65] | for the transition between levels A = 2p, A' = Is and taking into account that in this case the angular 
momentum of the photon can take values k = 1, 2, we receive four different terms contributing in Eq. ( |65] ). Unfortunately, in 
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the previous paper 1241] an error in the summation over projections of the Clebsch-Gordan coefficients did occur (see expression 
Eq. (41) in P24l]). Here we correct this mistake and give the proper expression for the probability. After angular integration and 
summation over projections we find 



where 



and 



h{L0) 



\h(L0') + /2(c^)r |J^(^) ^ ^^^JJ 



duj, 



1 





dridr2rlrle ''^ ^ gi{EA - U}\ri,r2) 



(66) 



(67) 



1 

71 



(68) 



respectively. Inserting again the representation Eq. f5% for the CGF with I = 1 leads to radial integrals that can be evaluated 
analytically. 

Substituting the integrals ( |67] | and ( |68] | into ( l66b and integrating over frequencies uj yields 



= ^ y" dH/j5piE2 ^ 1 98896 . iQ-5 (q,^)8 ^ 6.61197 • 10"^ s"^ {Z = 1) 



(69) 



with luq = i?2p — i^is- In Eq. i6% we indicated the Z-dependence of the W^^J^^ transition probability. Compared with the 
relativistic result in the "length" gauge (see 1 15i1 ) the relative discrepancy is about 0.1%. 

The calculation of the E1E2 two-photon decay with the nonrelativistic "velocity" form is more involved. Now the gauge 

constant should be chosen either K = — y/ -^ii for the form Eq. ( |45] | or X = for the form Eq. ( l44b . 
We choose K = together with the form Eq. (l44l i. The potential in this case reads 



fc + 1 



(2fc + 1)!! y fc(2fc + 1) 



fcri';K)^ + -Vfc(fc + i) 



d 



dr 



(70) 



The formula for dWf^f^ follows again from Eq. ( |65l ). Performing angular integrations and summations over projections as 
discussed in previous cases now yields 



with radial integrals of the type 



LO^ \h{L0) + hiuj't + + l2{L0)\' 



dhj 



(71) 



1 

71 



dridr2r^r2 e ^'^ 



1 - 



9i r2 



d 2i 
dri ri 



(72) 



and 



/2H = 



1 

71 



dridr2r\r\e ^ 



r2_ 
2 



d 3i 
dri ri 



g2{EA - w;ri,r2) 



(73) 



together with parameters v = Z/ yJ—2{E2p — uj) and v' = Zj 2(i?2p — ^'), respectively. The integrations over r\ and r2 
lead to a rather lengthy analytical expression containing various combinations of notations similar to those in Eqs. (|67] | and (|68] |. 
The numerical evaluation yields finally 



0)0 

W^a'^p iT dM^J'pT (^) = 3.6896 • lO"*' {aZf a.u. ^ 1.227 • 10-^ = 1), 



(74) 



12 



where ujo = 3/8 a.u. This result differs from that obtained from relativistic calculations ifTill . ifisll by about 0.5%. This 
discrepancy is larger than the relativistic corrections (of the relative order a^) and can be attributed to the pure numerical 
uncertainties in both calculations. Note, that unlike the case of the "length" form, the negative-energy contribution is no longer 
negligible when the "velocity" form is employed. Therefore, the result Eq. ( |74l i does not coincide with Eq. ( |69] l and represents 
only the positive-energy contribution to W^^^^ in the "velocity" form. Correspondingly, we compare this result to the positive- 
energy contribution calculated in lfT4|] . lITsll . The negative-energy contribution to W^^]^^ in the "velocity" form for low Z values 
was evaluated analytically in ifisll . 



ElMl TWO-PHOTON DECAY 

For the mixed ElMl two-photon transition probability the expression ( |64] | should be replaced by 

{A'\V^\Lo)\N)iN\V^\Lj')\A) 



E 



^ En - Ea + io 

{A'\V^\uj')\N){N\V^^{lo)\A) 

En- Ea+ lu' 
{A'\V^\uj')\N){N\V^^{uj)\A) 



En~ Ea+ lu' 
{A'\V^\lu)\N){N\V^\uj')\A) 



En — Ea + uj 



doj. (75) 



Here V^^{uj) = |tj3/2^yW^ V^^{uj) = ^^iquj^/^ {jim^. + sia/„), i^i/j is the spherical tensor of the rank 1 with the 

spherical component (i.e. the spherical function Yim^), Mo = § is the Bohr's magneton, jiM^ ™d sim„ are the spherical 
components of the total angular momentum and the spin operator (spherical tensors of rank 1) of the electron. This choice 
corresponds to the nonrelativistic "length" form for describing the emission of the electric photons. Since the potential for the 
magnetic photon includes total angular momentum and spin operator, coupled wave functions with the set of quantum numbers 
iV = {nlsjm] should be used, i.e. 

Cl^nlsj^ = E C'lZ sm. {t)Y^] {^Msnr. , (76) 
mirns 

where Rni (r) is the solution of the radial Schrodinger equation and Xsm^ = 1/2) is the spin function. The magnetic potentials 
in Eq. ( fTST i do not depend on radial variables. Thus, only the intermediate state with nl — uaIa or nl — ua'Ia' wiU contribute 
to the probability in Eq. (fTST l. After performing angular integrations and summations over all projections one arrives at the 
expressions 

= (77) 

and 

3/8 ^ ^2 3 



As the final result we obtain 

Wt^s' ^-(-Y -^^{aZf a.u. = 9.6769 • lO'^ s'^ [Z = 1). (79) 
^P^' TT \3j 655360^ ' ^ ' 

Again the Z-dependence of the W^^Jf ^ transition probability is indicated. Comparison with the result of a fully relativistic 
calculation now reveals a discrepancy of about 0.1%. 
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E1E2 AND ElMl TRANSITION PROBABILITIES FOR THE 2P HYDROGENIC STATE EXPRESSED THROUGH THE 

PHOTON MOMENTUM AND POLARIZATION. 



In this section we are going to derive expressions which allow the analysis of the dependence on the directions and polarization 
degrees of freedom of the emitted photons. For this purpose we turn back to the generic 5-matrix formulation employing the set 
of quantum numbers defined by polarization vector e and wave vector k. Relativistic units are used throughout this section. 

The S'-matrix element of the two-photon decay process A A' + 2'-f for the noninteracting electrons is represented by 



Sa^a = {-i)^e^ j dxidx2 {'^A'{xi)rA;{xi)S{xi,X2h''A;{x2)Mx2)) 

where xi — (ri, ^i) and X2 — (i'2, ^2) are 4-vectors. Also 

00 _ 



at \ ^ f J iuJi(ti-t-,) 4'n\'- 



(l-iO)+wi 



(80) 



(81) 



is the electron propagator where the sum runs over the Dirac spectrum for the electron in the field of the nucleus, (a;) is the 
electron Dirac wave function. En is the electron energy. 



^g(A)g»(kr-c^t) 



(82) 



is the wave function of the photon characterized by the momentum k and polarization vector (/^, A = 1, 2, 3, 4), a; = (r, i). 
For the real transverse photons 



K(x) 



ee 



2{kr— ojt) _ 



Ae,k(r)e" 



Inserting Eqs. (ISTTi-dSSTl in Eq. dSOl l, integrating over time and frequency variables and introducing the amplitude Ua' A as 



427) ^ _2i:i8{Ea' + J - Ea)U^a'A 



we obtain 



'-^A'A — 



27re^ ^ {aAl^^)A'n{aK',\.')nA 

^ En -Ea+UJ' 



(83) 



(84) 



(85) 



where e is the electron charge. 

Taking into account photon permutation symmetry we define the transition probability as 



dWA' 



2nS{EA - Ea' - UJ- uj') 



^A'A + ^A'A 



2 dk dk' 

(27r)3 (27r)3' 



(86) 



Taking dk = oJ^ dr\.\f^duo and integrating over oj yields 



^ En-EA+Uj' ^ 



^^ oj'{Ea-E'^-u') ^ 

mAm^ 
2 



1 



(27r)3 

{<^K'M')A'n{aA*)nA 



En - Ea+UJ 



2jA + 1 

dn]s^ni^i duj' . 



(87) 



Here the summation over n abbreviates the summation over the whole set of quantum numbers njnlnmn of intermediate states. 
The sums over the projections of the total angular momentum of the final state A' and the averaging over the projections of the 
total angular momentum of the initial state A in Eq. dSTl l are also included. 
In the Pauli approximation (see Eq. (|9])) the dSTT i expression can be rewritten as 



dW^„ 



P(27)_ 4,Uj'{Ea — Ea' — Uj') 



A'A 



E 



1 



E 



[/PJe,k)[/P^(e^kO ^ ^ £/Pje-,kO£/P^(e,k) 
E -n - Ea+uj' ^ E -n- Ea + UJ 



mAirij^ 
2 



•ijA + 1 

dr\.]ji\^' dio' 



(88) 
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where U^,j^{e, k) is defined by Eq. (|9]). 

Within the Pauli approximation we have to take into account the two terms of the exponent expansion in the expression for 
the emission operator 



A^{e, k) = ((e*p) + z(e*[k x s]))^-*-- « — [^,e*r] + -[H, (e*r)(kr)] + - (e*[k x ([r x p] + 2s)]) 



(89) 



The first term in this expression represents the electric dipole moment of the emitted photon, the second one is the electric 
quadrupole moment, and the last term in Eq. ( [89] ) represents the magnetic dipole moment of the emitted photon. 

For the beginning we evaluate the two-photon El and E2 decay rate. The probability of this process can be written as 



1 



(2-)^ + 1 

^ uT'\e,mT\e'y) + t/ir (e,k)t/„T^(e',k') 



En ~Ea+Oj' 

uT\e', k')L.„THe, k) + )(e', l.')uT\e, k) 



(90) 



Er^ + Ea' + UJ' 



the notations C/^;„ and Ujj^-^ for the dipole and quadrupole photons are used in accordance with decomposition Eq. 
The Hamiltonian H in Eq. ( |89] | acts on the eigenfunctions and, therefore, we can rewrite the expression Eq. ( l90l l in the form 

,'3,. ,3 I 



,TTr(ElE2)/ i\ 



E 



4(27r)3 ^ 2j^ + l 
(e*r)A'„((e'*r)(nk'r))„A ^ ((e*r)(nkr))A'„((e '*r))„A 



En-EA + uj' 



'E 



(e^*r)^,„((e*r)(nkr)), 

En~ Ea+i^ 



A . , ((e'*r)(nk'r))A'„((e*r))„^ 



En -Ea+UJ 



En -Ea+lu' 

2 

dnknk'dw', 



(91) 



where the relation ui = Ea — Ea' — uj' holds. 

For the summation over polarizations we use again Eq. ( fTSl ) and for integrating over the directions of the emitted photons we 
employ the formulas 



The formula 



/■ /■ 47r 

dnni = / dnriinkni = 0; / dnniUk = —Sik, 

An 

dnniTijnkni = —{SijSki + SikSji + 5a5jk)- 
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^ |ee'|2 = 1 + (ntrik')' = 2 - [nk x nkf 



(92) 



(93) 



is also useful for calculations. We should stress that in our paper 111511 an unfortunate misprint does exist: Eq. (51) should contain 
a vector product. 

For the evaluation of the quadrupole matrix elements in Eq. W\\ it is convenient to make use of the identity 



^imi(n)yi2m2(n) = Yl 



LM 



' (2^i + l)(2/2 + l) ^0 LM y 

47r(2L + 1) ^iio;20<-^/imi;2m2-^-f^Afin;. 



Therefore, 



(e*r)(nkr) = yiy^ i,,e;^(nk)p2r'>^2M(nr). 



P1P2M 



(94) 



(95) 
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In this case the integration over angles can be provided by the (see 113311 ') standard relations 



4n 



(96) 



f ^„V* f^W /„\V t^^ — I (^^^ 1)(2^3 + l) ^;iO Wimi 

J "■nYi^^^(n)ri^„i^(n)ri^m3[T^) - d 47r(2?i + 1) ^i20ho^i2m2i3m3- 
Finally, integrating Eq. ( |9T| l over angles and summing over all projections of the angular momenta, we receive 



(97) 



jTir(ElE2)/ /^ 4 



(271)35232 



j dnk y"dnk,^|eef [w'2[/i(w')+/2M]'+w2[/i(w) + /2(w')]']dw', (98) 



47r 47r 



where /i (w), I2 (w) are defined by Eqs. (|67] |. (|68] |. After the summation over polarizations and integrating over photon directions 
the previous result Eqs. (|66T l and (|69] l are recovered. 

The next step is the evaluation of ElMl two-photon transition 2p ^ Is + 7(M1) + 7(E1) with the set of quantum numbers 
of e and k. Formula Eq. ( |90] l can be cast into the form: 



,„.(E1M1). I. 4^Uj'{Ea- Ea' - i^') 

d^A'A 



^ 2m + 1 



E 



(99) 



E 



En + Ea' + OJ' 



and, according to Eq. 



dT/F\,^ ^(w ) = e 



4(27r)3 2jA + 



(e*r)^,„ e'*[nk' X (l + 2s)] 



K - -Ba + w' 



(e'*r)^,„ e*[nkx (l + 2s)] 

^ En~EA+U} 



e'*[nk' X (l + 2s)]) (e*r)„^(ek) 

A' 71 



T S 

e*[nk X (i + 2s)) (e'*r)„A 

/ A'n 

En-EA + uj' 
2 

dnkHk'dw', 



(100) 



where 1 = [r x p]. 

All the matrix elements can be easily evaluated with the aid of Eqs. (|96]l-(l97]i. Finally, we arrive again at Eqs. (|77]l-(|79]l: 



(ElMl) _ 4 



2^ (2 



TT V3 



12 



(3/8 - ((3/8 - u' f + duj', 



(101) 



3/8 

Ml ^ i / dW^jlY" = --^{aZf a.u. = 9.6769 • 10"^-^ {Z = 1). 
2 J ^P^" TT 10935^ ^ ^ ' 



ElMl 
2p 



(102) 



Thus, in this section we have evaluated again the two-photon transition probabilites, E1E2 and ElMl, for the process of 2p 
state decay in hydrogen atom. In contrast to |24], calculations were performed based on the representation dealing with the 
other set of quantum numbers. The analytic results obtained are in a good agreement with corresponding relativistic values 
(discrepancy is not more then 0.1%). Recently the two-photon emission probabilities was evaluated by fully numerical methods 
in |44]. The corresponding results are in perfect agreement. The method applied in this section simplifies the determination of 
the dependence on directions of the emitted photons in the two-photon transition probability when the external field is present. 
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TWO-PHOTON DECAY OF THE 2S AND 2P HYDROGENIC STATES IN AN EXTERNAL ELECTRIC FIELD 



In this part of our work we evaluate the two-photon 2s- Is transition probability for the hydrogenic atoms in the presence of 
an external electric field. Similar as in section 11 of this paper we take into account only the mixing of the 2s and 2p states. Such 
level mixing leads to additional E1E2 and ElMl two-photon decays besides the dominant "pure" ElEl two-photon transition. 
As it was shown in previous sections the E1E2 and ElMl two-photon transition probabilities are about (aZ)'^ times smaller 
then ElEl transition probability. However, a consideration must be given to the fast increasing accuracy of the spectroscopical 
experiments. Therefore, E1E2 and ElMl two-photon transition rates should be taken into account like a correction to ElEl 
decay. Moreover, we will show that in the presence of the external electric field terms linear in the field will add to ElEl 
transition probability; These are the interference terms. Furthermore, most if not all spectroscopical experiments (see most 
accurate experiments 10, IH) involve an external electric field so the analysis of the two-photon decay transitions in an external 
electric field is indeed required. 

In order to evaluate the two-photon 2s ^ Is + 27 transition probability in an external electric field we should turn to Eq. (ISST i 
and use the decomposition Eq. ( [89] l. The mixing of the 2s and 2p states is described by Eq. (fT6] l. According to Eq. (|20] |. the 
two-photon transition probability in an external weak electric field can be introduced in the form: 

' „ . , 4 ^''^ V- 1 1^ (ls/.'|AP(e,k)|n)(n|AP(e^kO|2sM") 

-3„i:i:(-i)w;£.... '"'''|-^^'t'^'i"'f!r''''"^"" + 



E 



En - Ea+Oj' 

(lsAi'|ylP(e', k')|n) {n\A^{e, k)|2s/x") 



\-\-( nr^-^ (lsA.1^P(e^kO|n)(n|AP(e,k)|2pM) 



En - Ea + 

2 



n fiq 



where (e, k) is defined by Eq. 

The decomposition Eq. ( [89l ) of the (e, k) operator shows that the first and third terms in Eq. ( I103l l correspond to the ElEl 
two-photon decay rate of the 2s electron level. Other terms in Eq. ( 11031 ) represent the admixed E1E2 and ElMl amplitudes of 
the 2p state two-photon decay probability. Application of the expansion Eq. (|89] l to the second term in Eq. ( |103l l leads to the 
expression 

^ {ls^i'\AP{e,k)\n){n\APie',k')\2p^i) . „^ (ls/i'|er|n) (n| l(e^r)(nk^r)|2pA.) 

^ En-EA+Lo' ^ En-EA+LO' 

(ls^'|e*r|n)(7i|i ([e' x nk'](i + 2s)) \2p^I) 
+ujuj'22 ^ ^ ^ — : — ; (104) 



■ / 2 \ ^ 

—luj ui y 



En-EA+ to' 
(ls/i'|i(er)(nkr)|n)(n|e'r|2p/i) 



En- Ea+ uj' 
{ls^l'\\ ( [e X nk](i + 2s)) |n)(n|e'*r'|2p^) 

En-EA+ LO' ■ 



In order to perform the summation over intermediate states in this expression we employ the Coulomb Green function method, 
which was described in the third section of this paper We should note that in terms involving a magnetic dipole photon the 
Coulomb Green function does not occure due to the orthogonality of the radial functions. 
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Summing over aU angular momenta projections in the expression ( 11031 ) we obtain 



UJLU 



2 '2 |„„/|^ r7-2s Is/ N I 7-2sls/ / 



3-5 A2 



2^ ( I) ''(^) + (e [eD x [e' x nk']]) + uj' (e' [eD x [e x nk]])] (ee')' 



3 J A2 

^^2^,2 2V3^ |eef [/f + 7?^ ^^(c.')] [C./2'" '^(^')(eDnk) + c.^^^ ^^(^)(eDnkO] 



+ 



2-3 
'52A2 
2-3 



52A2 



(ee')(DnkO^'/iP''(w') + (ee')(Dnk)t^/iP ''(w) 
(ee')(Dnk)w/2P''(t^') + (ee')(DnkO'^'/2'' ''(w) 



4-3 

'52A2 



(ee')(Dnk'V/i (^') + (ee')(DnkV/rP ''(c^) {ee')i-Dn^)ujI^P '\uj') + {ee'){T>ni,,)uj' if" (lj) 

■ ^[c^([Dxe][nk' xe'])+c^'([Dxe'][nkXe])]' 



215 



37 A2 



(De) [nk' X e'] + uj' (De') [nk x e]]' 



dnkdnkfduj' , (105) 



where ^^{^) is defined by Eq. ( |60l ) and Eq. ( |6T| i for the 2s-ls transition, and I^^ ^^{^), ^^{^) ^re defined by Eqs. ( |67| i. 
(l68T l. correspondingly. 

The summation over polarizations is provided by Eqs. (fTsT l and ( |93] l. The radial integration can be analytically performed as 
well. Finally, the integration over ui' is carried out with MATHEMATICA code. The result is 



(27) 



- 0.00131822(.Z)«-^:^1^^^^^'^' 



dnkdn. .3 ^2 4Dnk + Dnk'](l + (nknkO^)(oZ)^ 

- "■~^' ^[Dnk + Dn.] (1 + (nkn.)) (aZy + 0.0017509l4?(.^)^ 



A2 



(106) 



The first term in this expression is the differential ElEl 2s-ls transition probability. The second and third terms represent 
interference terms for the mixed 2s and 2p two-photon transition probabilities (second term corresponds to the E1E1E1E2 
transition probability and third one represents ElEl-ElMl transition probability). Last term is related to the sum of the E1E2 
and ElMl two-photon probabilities, respectively. 

Rewriting Eq. (1106b in a form similar to Eq. (l22l l yields 



dW-, 



(27) 



2s Is 



dnkdnk' 



= 1^0 [1 ± MD) [nonk + nnnk'] (1 + (nknk')') ± ^2(0) [nnnk + nnnk'] (1 + (nknkO)] , (107) 



where Wo = W^^'^^ + W^l'^'> e'^ D'^ / , W^l'^^ is the sum of the E1E2 and ElMl transition probabilities. 
Functions (3i{D) and /32{D) are defined by 



P2{D) 



0.000230135(aZ)^ [elLT 



2p 



W^n7r3 A2 ' 

0.0000340919(aZ)^ lelDFap 



(108) 



1^071-3 

The maximum value for (3i (or (32) is achieved at the field strength 

w27A 



A2 ■ 



±0.000018 a.u. w ±57kV/cm, 



(109) 



where the (— ) and (+) signs in Eq. dllOl l correspond to the H and H atoms, respectively. 
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Then the corresponding maximum value of the transition rate dW-^'l'' is obtained via 



dW-, 



(27) 



2s Is 



M^o(i^max) 



X [1 ± 0.00024397 [nonk + nnnk'] (1 + (nknk')^) ± 0.00003614 [nnm, + nnnk'] (1 + (nknk'))] ■ (HO) 

Integration over the directions of the emitted photons nk and nk' leads to the value Eq. {6% for the ElEl two-photon decay 
rate of the 2s-state, because the interference terms give zero result (see Eq. (|92||). But for the differential transition probability 
Eq. (I107l l the interference terms clearly demonstrate the linear dependence on the external electric field D. Thus the total 
two-photon W^'J^ transition probability integrated over photon directions is: 



WoiD„ 



Di 



A2 



3.98116 • IQ-^*^ a.u. w 16.4585 s" 



(111) 



i.e. twice as large as the zero field value Eq. ( fSSl ). 

In principle, the dependece on the external electric field in the transition probability W^^^^^ Eq- (111 II) can be consid- 
ered as a correction which does not vanish after integration over the photons emission directions. If we return to the ra- 



diative correction considered in |12], then it is easy to see that the radiative correction (Eq. (36) in ll2ll ) i5r2s/r2s = 
2.020536f (aZ)2ln [{aZ)-"^] ^ -2.4594 • 10"^ corresponds to the magnitude of the field |A| ~ D^^^y/sf^IJ^s ~ 
2.8 • 10"^ a.u. sa 90 V/cm. Such fields are often used in the spectroscopic experiments, therefore, this effect also should be 
included in this context. 

The linear over field corrections Pi{D) and ^2(0) in Eq. ( llOTI i reach the magnitude of the radiative correction at the fields 
approximately of the same order Unlike the correction, discussed earlier in the section "2s decay rate for hydrogen and anti- 
hydrogen atoms in external electric fields", this is the correction directly to the same process, as radiative correction [ 12]. 

As it was mentioned above the formal T-noninvariance of the factor nnnk in Eq. (ISTT l and in Eq. (1107b (nk and hd are T-odd 
and T-even vectors, respectively) is compensated by the dependence on r2p; This is the imitation of T-noninvariance in unstable 
systems, as predicted by Zeldovich [39]. 

The relative difference for the decay rates in H and H atoms at the maximum value Dmax equals to: 



2/3i(A, 



VFoPmax; 

^2^ 



)dnkdnk' Tyo(^max)rfnkdnk' 
,)(nDnk + nDnk')(l + (nkUk')^) + 2/32(£'max)(nDnk + nDnk')(l + nkUk') 
(nnnk + nDnk')(0.000280111 + 0.00024397(nknk') + 0.0000361414nknk')- 



(112) 



This ratio is close to 0.028% and represents a tiny effect reflecting the difference between matter and anti-matter even at maxi- 
mum field strength Umax- 

For completeness the 2p ^ Is + 27 two-photon transition probability should be considered as well. It can be evaluated 
similarly to 2s ^ Is + 27 two-photon decay rate with the use of the wave function 



|2p/i") - |2p/i") -7?^(2sAi"|eDr|2p/.)|2s/i" 



(113) 



In this case the two-photon transition without an external electric field will be provided by the sum of the E1E2 and ElMl 
decays, and the interference terms will be the same as in Eq. (11051 ). 
The result can be presented in the form 



2p Is 



A2 



0.0000340919 r2p 

^3 A2" 



eDnk + eDnk'] (1 + nknk') {aZf 



dnkdnk', (114) 



2£i^ = Wo [1 T l3i{D) [nonk + nnnk'] (1 + (nknk')') T HD) [nonk + nonk'] (1 + nknk')] , (115) 
ankunk' 



where Wq = W^^lf^'^ + W^^"^^'"' + 9e^D^W^^l^"' / and the functions Pi{D), (32(0) are defined again by Eq. (fToSl) 



^(ElMl) 
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Then the maximum of the (3i (or P2) is achieved at 

,„ , A 



3|e|w27 



7.1 • lO^^^a.u. « 0.23V/cm. (116) 



The corresponding maximum value of dW^ll IS 



^"'^ =l^o(Anax) 



X [1 T 0.00048613 [nont + nonk'] (1 + (nkUk')^) T 0.000720147 [nnnk + nonk'] (1 + nkUk')] ■ (117) 

After integration over nk and nk' in Eq. (II 14l i the term quadratic in the external electric field for the correction to the two-photon 
transition probability 2p-ls still remains: 

Finally, the (+) and (— ) signs in Eq. (II 171 ) correspond to the H and H atoms, respectively. The relative difference for the decay 
rates in H and H atoms at the maximum value Dmax equals to: 

2p Is /jjN 2p Is (H) = (119) 

W^o(^max)c^nkdnk' Wo(^max)c^nkdnk' 
= 2/?i(Anax)(nDnk + nDnk')(l + (nkHk')^) + 2/32(Anax)(nDnk + nDnk')(l + nkHk') 
= (nnnk) + nDnk')(0.0111629 + 0.0097226(nknk')^ + 0.00144029nknk'). 

This ratio turns out to be close to 1%. However, any direct observation of this difference should be difficult due to the huge 
background from the one-photon transition 2p — > Is + 7. 



THREE-PHOTON DECAY RATE FOR THE 2p STATE OF HYDROGEN-LIKE LIGHT ATOMIC SYSTEMS. 

In this last section we present the calculation of the ElElEl transition probability for the 2p hydrogenic state decay. Parametric 
estimate can be easily obtained in usual way and it is a{aZ)^ a.u. Therefore, one can expect a numerical result for the ElElEl 
decay rate which is comparable with values obtained for the E1E2 and ElMl transition. 

According to the Feynman rules the 5'-matrix element for the n-photons emission process with the transition from the state A 
to A' (A — > 717 + A') can be written as (in ru.): 

'^A/A = J (i'A'ixi)A{xi)S{xi,X2)A{x2).-.A{Xn-l)S{Xn-l,Xn)tpAiXn)^ dxi...dXn (120) 

where A is the emission operator, ipA, iPa' ^re Dirac wave function, S{xn-i, Xn) etc. denotes the bound electron propagator, 
and e is the electron charge. 

The emission operator is given by 

A{x) ^ , &(^)e-^('''-'^*\ (121) 
V (jj 

where e^'^' is the 4-vector of the photon polarization (e'^-' = 67, e is the polarization vector, 7 are the Dirac matrices) and k is 
the wave vector In what follows we employ the nonrelativistic approximation, replacing the Dirac solutions for atomic electron 
orbitals by the corresponding Schrodinger ones and omitting, where it is justified, the contributions of the negative energy states 
to the exact electron propagators in Eq. ( 11201 ). 

Integrating over time variables in Eq. ( 1120b results in 

Sy^ ^ -2Tn{-ie) — S{Ea' + nuj - Ea) j- — 7, (122) 

[Es„_i - Ea - (n- 1)lu\ ... [Es^ ~ Ea~uj\ 



UJ 



Si ,. .. 



20 



where the summation over Si extends over all intermediate Schrodinger states with positive energy, Es are the Schrodinger 
energies for an atomic electron. 

In the nonrelativistic approximation we can expand the exponents e*"^"" in Eq. ( 11221 ). leaving only the first term and replace 
the matrix elements {e)ss' with Dirac wave functions by the matrix elements (ep)^^/, with Schrodinger wave functions, where 



p is the electron momentum operator Then, using a well known quantum mechanical relation I [r, J 
amplitude of the three-photon emission process can be written in a form: 



i{p)ss', the 



E 

.Sl,S2 

+ E 

Sl,S2 

E 



Ua'\ = -(-je)3(27r)3/Vwit^2t^3 S^Ea' + loi + lo^ + UJ3 - Ea) x 
(e|ri)^,,^ (621-2),,,, (631-3),,^ , (etri)^,^^ (egrg),,,, (621-2),,^ 

{Es,_ -EA+0J1+ LU3){E,, -Ea+ UJ2) 

(Es^ --EA+LU2 + LU3)iE,, -Ea+uji) 
{e*3r3)A's2 (621-2),,,, (e^ri),^^ 



{E,, -Ea+uji+ uj2){Es, -Ea + UJ3) 

(e^ra)^,,, (e^ri)^^^^ (egrg),,^ 
{E,, -Ea+uji+ uj2){.Es, -'Ea + UJ3) 

(631-3)^'.. (e^n),,,, (e^i-2),^^ 



E 

Sl ,S2 

E 



(123) 



E 



^ ^ {Es, ~Ea+uji+ LU3)iEs^ -Ea+ UJ2) ^ {Es, -Ea+uj2+ W3)(S„ -Ea+ uji) 

Here wi, L02, ^^3 are the frequencies (energies) of the emitted photon. 
The probability of the three-photon emission process is 



(3) _ 



A' A 



2-K 



'^A'A 



d{EA' + LUi + u}2 + UJ3 - Ea) 77^^^^ . 

[Ztt 



(124) 



The total probability can be received from Eq. ( I124l i by the summation over photon polarizations ei, 62, 63 and integration over 
all the photon-emission directions ki, k2, ka and frequencies wi, 0^2, ^3. For the summation over polarizations and integration 
over photon directions it is convenient to use the relations Eqs. ( fTSl l. (|92| i. (|93] l. 

After averaging over angular momentum projections of the initial state and summing over final ones, the probability of the 
three-photon emission process results as 



9o4, ,3, ,3, ,3 



du}idu}2duj2 



5{Ea> + + L02 + UJ3 - Ea) X 



(125) 



E E (-l)'^^'^^'^^i'A(9l,'?2,g3)f/i^!4*(-9l,-92,-g3), 



f 3) 

where the expression for the amplitude U\,j^{qi, (721 93) in spherical representation is given by 



f^i^!4('?l,92,93) 



E 

.S1S2 

+E 

S1S2 

+E 



(('^i)'?i)a's2 ((''2)g2)„,i ((1-3)93),!^ 

{Es2 -Ea+LUi+ UJ2){Es, ~Ea+ LO3) 

((r2)g2)A'^2 ((''l)gi)s2Si ((1-3)93),,^ 

{Es2 ~Ea+lui+ oj2){Es, -Ea+ L03) 

((i-3)93)a's2 (('■i)9i)«2«i (('■2)92),iA , 



{Es, -Ea+lui+ uj3){Es^ -Ea+ L02) 



^E 

S1S2 

^E 

S1S2 

E 



ii^l)qi)A's2 ((l"3)93),„i ((1-2) 



Q2JsiA 



{Es2 -Ea+oji+ u;3){Es, - Ea + UJ2) 
((r2) 

<12)a'S2 ((l-3)93),2si ((l"l 

{Es2 -Ea+U)2+ i^3){Es^ -Ea+ COi) 

((i-3)93)a'.2 ((^2),2),,,i ((l-l)9l),iA 



{Es, -EA+OJ2+ ^3){Es, -Ea+ LOi) 



(126) 



In order to calculate the transition probabilities for the process 2p — > 37(E1) + Is in the hydrogen atom the nonrelativistic 
Coulomb Green function method is employed. Inserting the Green function in a form Eq. (ISOl l in (I126I I and representing the 

vector component (r), like = \/ ^Yiq we can perform the angular integration, which gives 



UA'Aill'l-i^Qs) 



21a + 1 
21a' 



1 2^ 2^ ^Wh^'-'WhO^WlAGX^l^'-'^^^'^'^'^' y~'lqil2nii,^lq2hrni^ 



a 



I133 Iaitia 



hmi hrrii 



^lq2 hmi, Iqi hmi^ I93 lArriA 



.hmi^ 1 I p (,, ,, \ [/^'a'^a' f^<'2m.i, „lirni^ ^Ij^im^, ^hmi, „iinili 

+ ^hh l'^3, j [L-iqj i,mi, '-'Iqa hmi^ Iq, UmA ^Iqs /2m,, "^Igi /im,^ "^192 

-rri^l,yiy5,iye) ^'-^iq, i,m,^^iq^ i^rn,-^^lqi UrriA '-'I133 bm,, ^1^2 iim,^ '"'Igi Z^mAj J 



(127) 
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Here c!^™''* , ^ is the Clebsch-Gordan coefficient, and 



oo oc oo 



Fiuhi^t^'^j) ^ J J J dridr2dr3r^r2r^Rn^,i^,{ri)gi2{'^i;ri,r2)gii{i^j;r2,r3)R„^^i^{r3), (128) 





lyi = Z/^-2{Ea-uji-u;2), 1^2 = Z/^-2{Ea- l^s), = Z/^-2{Ea - - oJz) (129) 
z/4 = Z/^J-2{Ea-uj2), = Z/v/-2(Sa-c<J2-^3), = ^/v/-2(£;a - Wl) 

corresponding to the 6 different terms in Eq. ( 1126b . 

For the 2p — > 37(E1) + Is process, Rua.Xa = ^2p, R-n^'iA' = ^i*" ~ ^' ~ ^- Performing the summation over all 
angular momentum projections in Eq. ( 11251 ) we arrive at the expression 

dwgis = e'iuMcol^^^^^i^piE,^ + + + C.3 - i?2p) [l5Foi(^i, 

+ 15Foi(z/3, 1/4)^ + 15i^oi(l^5, I'e)^ + 20Foi(l^5, i^6)i^2l(i^l, 1^2) + 12^^21(1^1, ^^2)^ + 20i^oi(l^5, i^6)i^2l(i^3, 1^4) 
+4F2i(z^l,J/2)^2l(i^3,l'4) + 12F2i(i^3,J/4)^ + 4i^21 (i^l , 1^2)^21 (i^5 , I'e) + 4^21 (1/3 , Z^4)-F21 (i^5 , I'o) 

+ 12F2i(l/5,J^6)' + 10i^Ol(j^3,J^4){Foi(!^5,J^6) + 2i^^2l(i^l,J^2) + 2F2i(l/5,J^6)} (130) 

+10^bi(i^i, t^2){i^oi(i^3, ^■4) + -Foi(i^5, J^e) + 2i^2i(i^3, i^i) + 2^21(1/5, i^e)}] ■ 

After performing the radial integrations we should integrate over frequencies. The (5-function in Eq. (I130l l annihilates one of the 
integrations. It is convenient to perform the other integrations after the transformation of the variables: ui = jAE{l — x){l+y), 
LU2 = |Ai5(l - x)(l - y) and lo^ = ^'^E{1 + x). Here AE = E2p - Eis = 3/8 in atomic units. The final result will be 
written in atomic units; 

3/8 3/8-U3 3/8-^3-1^2 3/8 3/8 3/8 

W^l\^ = Jdu;^ J du;2 J du;,dW^X ^hj J J '^^^^^'^ ^ 
00 000 

2 1 1 

= ^ (i) / / = 0-263466(1) . 10-^a(aZ)«a.u. = 6.391(1) • 10"^ (131) 

-1 -1 

For Z = 1 the result is VK2p\s — 6.391(1) • 10^^ s^^. The smallness of this result even compared to the two-photon 2p-ls 
decays means that the 3-photon transitions cannot play any significant role for the astrophysical purposes, i.e. for the electron 
recombination history in the early Universe. For comparison we also provide the one-, two- and three-photon transition rates for 
transitions from 2p, 2s levels (in ru. and s^^, where m is the electron mass): 
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TABLE I: Electric (El) and magnetic (Ml) multi-photon transition rates between 2s, 2p and Is levels are presented for hydrogen-like ions in 
relativistic units (r.u.) and for hydrogen in units s~^, respectively, m and Z denote the electron mass and nuclear charge number. 



transition (type) transition rates 



2p 


Is 


(El) 


3.902 


10~^ mQ(aZ)'' 


r.u. 


6.25 


■ 10** s- 


1 


2s 


Is 


(ElEl) 


1.317 




r.u. 




8.229 s- 


1 


2p^ 


Is 


(ElElEl) 


2.635 


IQ-^ ma^{aZf 


r.u. 


6.39 • 


10-*s- 


1 


2s ^ 


Is 


(Ml) 




^ma{aZy° 


r.u. 


2.5 ■ 


10"^- 


1 


2p^ 


Is 


(ElMl) 


2.911 


10-5 ma2(Q^)8 


r.u. 


9.68 • 


10"^- 


1 


2p^ 


Is 


(E1E2) 


1.989 


10-^ ma^iaZf 


r.u. 


6.612 • 


lO-S" 


1 



CONLUSIONS 

In this paper analytical results for 2s, 2p levels decays for the hydrogen-like atomic systems with one-, two- and three- 
photon emission have been presented. All calculations were performed within Pauli approximation utilizing the Coulomb Green 
function. The emission processes were also calculated in the presence of an external electric field. For the two-photon decays in 
the absence of an external electric field the obtained results are in good agreement with those of other calculations. 

The two-photon decays E1E2 and ElMl were evaluated with different sets of quantum numbers (representations) for the 
emitted photon, namely, parity and momentum or polarization vector e and wave vector k. Moreover, we have employed 
different forms in combination with different gauges. The results do not differ in magnitude by more then 0.1% from fully 
relativistic values which were obtained earlier. Recently a paper |44] did arrive where, in particular, 2p-ls ElMl and E1E2 



transition rates were evaluated for the H-like ions in the wide range of Z values. For Z = 1 the results of 114411 are in agreement 
with our ones. 

We have evaluated also the two-photon decay probabilities E1E2 and ElMl with the set of quantum numbers e, k to investi- 
gate the probability dependence on directions of the photons emission. It allowed us to obtain the two-photon emission processes 
in the presence of the external electric field. In particular, we have demonstrated that interference terms in the ElEl and E1E2, 
ElMl transitions appear, which depend linearly on the external electric field. 

The important result of our calculations is the prediction of a characteristic difference in transition probabilities (spectra) 
between H and H atoms in the presence an external electric field, caused by the terms, linear in the electric field. This effect 
was not yet discussed in literature. The observation of this effect would also allow for drawbacks on our understanding of 
fundamental symmetries in nature, i.e. the CPT-symmetry: Any deviation from this result would provide a hint for CPT violation 
from a low-energy physics scenario, i.e. atoms in an external electric field. 

Finally, we have compared the radiative correction evaluated in |12] with the electric field correction and determined the 
magnitude of the electric field strength, when both corrections become of the same order 
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